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We present the experimental investigation of the non-Gaussian nature of some mixtures 
of Fock states by reconstructing their Wigner function and exploiting two recently in- 
troduced measures of non-Gaussianity. In particular, we demonstrate the consistency 
I between the different approaches and the monotonicity of the two measures for states 

belonging to the class of phase randomized coherent states. Moreover, we prove that the 
I exact behavior of one measure with respect to the other depends on the states under 

• • investigation and devise possible criteria to discriminate which measure is more useful 

^ for the characterization of the states in realistic applications. 

• 

Keywords: Quantum optics; Photon statistics; Quantum state engineering and measure- 
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1. Introduction 

Gaussian states, i.e. states with Gaussian Wigner functions, are the key ingredi- 
ent of many continuous variable (CV) Quantum Information protocols GE How- 
ever, in order to achieve some relevant tasks, non-Gaussianity (nonG) in the form 
of non-Gaussian states (states endowed with a non-Gaussian Wigner function) or 
non-Gaussian operations is either required or desirable. For instance, it has been 
recently demonstrated that nonG can be used to improve teleportation, cloning and 
storage; in addition non- Gaussian operations are interesting for the realization of 
entanglement distillation and noiseless amplification P 

Several implementations of non-Gaussian states have been reported so far, in 
particular from squeezed light, 33 close-to-threshold parametric oscillators^ and in 
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superconducting circuits!^ Such states have been mainly achieved in the low en- 
ergy regimdSEES by using single-photon detectors, visible light photon counters 11 
and time-multiplexed photo-resolving detectors.^ More recently, we extended the 
investigation to the mesoscopic regime by exploiting the linear response of hybrid 
photodetectorJUl 14 ! 15 ! and Si-photomultipliersP^ 

Being it recognized as a resource for CV Quantum Information, the need of 
quantifying the non- Gaussian character of states and operations naturally arises and 
different non-Gaussianity measures have been proposed. 17 Indeed, not all of them 
are characterized by an operational meaning. Moreover, in the realistic situations, 
non-Gaussianity measures based only on quantities that can be experimentally ac- 
cessed are desirable. To this aim, here we present an experimental work in which 
we compare the two measures introduced in Refs. [18] and [19] by testing them on 
phase-randomized coherent states or phase- averaged coherent states (PHAVs), a 
class of states exploited in communication channels and in decoy-state-based quan- 
tum key distribution protocols P2 The two measures, being in perfect agreement 
with each other, can be considered a useful tool to quantify the non-Gaussian na- 
ture exhibited by the Wigner functions and testified by the experimental results. 
Moreover, the sensitivity of such non- Gaussian measures to small changes in the 
mean number of photons suggests their possible exploitation to test decoy states 
based on the class of PHAVs. 

2. Quantifying nonG amount 

A single-mode PHAV gp is a classical state obtainable by randomizing the phase (j) 
of a coherent state (3 = \/3 \ e 1 ^ . It is characterized by a density matrix diagonal 
in the photon-number basis and by a Poissonian photon statistics. The state is 
obviously phase-insensitive and its Wigner function reads 

Wphav(«; p) = ^ e-l^l 2 = - J (4M|/3|) exp[-2(|a| 2 + \/3\ 2 )] , (1) 

where Iq(z) is the modified Bessel function. This function, being it endowed with 
a dip in the origin of the phase space, is clearly non-Gaussian. 

Another state exhibiting a nonG nature, which can be useful for application to 
passive decoy state quantum key distribution,^ can be obtained from the interfer- 
ence of two PHAVs gp and g„ (2-PHAV hereafter) at a beam splitter (BS) with 
transmissivity t. The states outgoing the BS are still diagonal in the photon number 
basis and the transmitted one is characterized by the following Wigner function: 

W 2 -PHAv(a; P,M= [ ^ WpuAv(a - pVl^t; M) (2) 
Jo 27r 

where ft = \(3\ e 1 ^ and the function in the integral is given by Eq. |l]). Obviously, 
the Wigner function of the reflected mode can be obtained by replacing t with the 
reflectivity (1—t). In both cases, the quantification of nonG amount can be achieved 
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by considering two competitive measures recently introduced. One of them is based 
on the Hilbert-Schmidt distance from a Gaussian reference state, namely: 

5M = T^r = « ' (3) 

where ji[g) is the purity of the state g : a is a reference Gaussian state with the 
same covariance matrix as the state gp under investigation and K,[Qp,a] =Tr[gpa] 
denotes the overlap between gp and crP^l 

The second measure we address is given by: 

&B[Qp] = S(v)-S(Qp), (4) 

where S(g) = —Tr[g\ng] is the von Neumann entropy of the state g^ As both 
PHAV and 2-PHAV are diagonal in the photon- number basis, their reference state 
is a thermal equilibrium state, with the same mean number of photons N = |/3| 2 . 
Moreover, the two measures result based only on quantities that can be experimen- 
tally accessed by direct detection, since both Eqs. ([3| and (4) can be expressed in 
terms of photon-number distributions! 2 ^ In particular, Eq. 

n 1~ n ) 



3|) reduces to 

(5) 



En^ 2 

where r n = N n / (1 + 7V) n+1 is the photon-number distribution of a single-mode 
thermal state having N mean number of photons and p n is the statistics of gp. In 
particular, PHAV is described by a Poissonian distribution 

pP HAV =exp(-|^| 2 )|/?| 2 "/n!, 

whereas 2-PHAV is characterized by a non-trivial 2-peaks distribution! 23 ! 24 ! 

„2-phav_^ P -A^fn\ (-i) fc /i?\ fc r(l/2 + fc/2)r(l/2) 

n\ f^ \k) 2tt \AJ T(l + fc/2) {> 

x X F 2 [{1/2 + fc/2} , {1/2, 1 + k/2} , B 2 /4] , 

in which A = \/3\ 2 + \/3\ 2 , B = 2|/3||/3| and \F2(a,b,z) is the generalized hypergeo- 
metric function. 

On the other hand, Eq. Q becomes 

Ss[Qp] = (N+l) ln(N + 1) - NlnN + \np n . (7) 



3. Experimental results and discussion 

We generated the class of PHAVs by exploiting the second harmonics (@ 523 nm, 
5-ps pulses) of a mode-locked Nd:YLF laser amplified at 500 Hz (High-Q Laser 
Production). According to the experimental setup sketched in Fig. [I] we obtained 
the single PHAV by sending the light pulses to a mirror mounted on a piezo-electric 
movement. Its displacement, which was controlled by a function generator, was op- 
erated at a frequency of 100 Hz and covered a 12 /im span. Moreover, we produced 
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Fig. 1. Experimental setup. Fj: variable neutral density filter; BS: 50/50 beam splitter; Pz: 
piezoelectric movement; MF: multimode fiber (600 jim core). 



the 2-PHAV from the interference of two single PHAVs at a BS. A continuous vari- 
able density filter Fi allowed us to change the total energy of the states, whereas 
a second filter F2, inserted in the path of one of the two PHAVs, was used to 
change the balancing between the two components of the 2-PHAV. As the states to 
be characterized can be fully described by their photon-number distributions, we 
implemented a direct detection scheme involving a photon-counting detector, i.e. 
a hybrid photodetector (HPD, R10467U-40, maximum quantum efficiency ~0.5 at 
500 nm, Hamamatsu). This detector is characterized not only by a partial photon- 
counting capability, but also by a linear response up to 100 photons. Thanks to 
its properties, the HPD can actually operate in the mesoscopic domain, where the 
states are robust with respect to losses. The output of the detector was amplified 
(preamplifier A250 plus amplifier A275, Amptek), synchronously integrated (SGI, 
SR250, Stanford) and digitized (AT-MIO-16E-1, National Instruments). The gain of 
the detection apparatus was obtained in a self-consistent way without any a priori 
calibration! 25 ! 26 ! This method allows us to reconstruct the detected photons distri- 
butions of the states, which represent the basic element to retrieve their Wigner 
function. Such a goal can be achieved by mixing at a BS the state to be charac- 
terized with a coherent probe field whose amplitude and phase are continuously 
varied EH 28 1 29 1 In this case, as both the states are phase-insensitive, we actually 
reconstruct only a section of their Wigner distribution for fixed phase. 

In Fig. [2] we plot the experimental data (blue dots) of a single PHAV (left 
panel) and of a balanced 2-PHAV (right panel), endowed with nearly the same 
mean number of detected photons, Mr = 1.97 and Mr = 1.94, respectively. In 
each panel we also show the 3D-theoretical expectations (mesh) for the PHAV and 
2-PHAV, respectively!^ 

WpRAv(Vta) = W P KAv(Vta)e-^(\«W\)^ (8) 

£ being the overall (spatial and temporal) overlap between the probe and the PHAV, 
and: 



WVphav(vW) = W 2 -PHAv(v / Cp«)e- [v ' I3??|a|+v ^ S( ^ l+ ^ l)1 , (9) 
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Fig. 2. Experimental reconstruction of a section of the phase-insensitive Wigner function of a 
PHAV (left panel), with \(3\ 2 = 1.97 and f = 0.999 and of a 2-PHAV (right panel), with |/3| 2 = 1.03, 
|/3| 2 = 0.91, £p = 0.95 and £s — 1- Blue dots: experimental data; mesh: theoretical surface. 

where £p describes the overall overlap between the probe and the 2-PHAV and £s 
the overall overlap between the two components of the 2-PHAV. In Eqs. ([8| and 
([9]), |/3 1 2 and |/5 1 2 are now the mean numbers of photons we measured (see Fig. [2]), 
thus including the quantum efficiency. In fact, it is worth noting that for classical 
states the functional form of the Wigner function is preserved also in the presence 
of losses and its expression, given in terms of detected photons, reads W(a) = 
f Xlm=o( — ^) m Pm l ai where Vma, represent the detected-photon distributions of the 
state to be measured displaced by the probe fieldpS As testified by the very high 
values of the overlaps £, £p and £s reported in the caption of Fig. [2j we actually 
achieved a very good superposition in aligning the system. From a direct comparison 
between the two panels it emerges that the states under investigation are non- 
Gaussian. In fact, the Wigner function of a single PHAV is characterized by a dip 
in the origin of the phase space, whereas that of a 2-PHAV with almost the same 
mean value exhibits a peak in the origin followed by a "shoulder". Moreover, it 
is worth noting that the measurements were actually performed in the mesoscopic 
photon-number domain, as the reconstruction of the Wigner functions was achieved 
by displacing either the PHAV or the 2-PHAV with a coherent field whose intensity 
was changed from zero up to four times the mean value of the states themselves. 

To quantify the nonG amount, we considered the measures introduced in Sec. [2] 
As from the experimental point of view we do not have access to photons, we 
calculated similar expressions, ca and 6b, for detected photons, which represent 
lower bounds to nonG EES In particular we found ca = 0.207 ± 0.004 and €b = 
0.156 ± 0.020 for the single PHAV, whereas we obtained e A = 0.036 ± 0.005 and 
6b = 0.012 ± 0.025 for the 2-PHAV. The consistency between the two measures, 
together with the fact that measure €a can be directly expressed in terms of Wigner 
functions^, demonstrate that a Wigner function exhibiting a dip in the origin of 
the phase space is more non-Gaussian than one characterized by a peak in the 
origin followed by a "shoulder". Moreover, the results prove that combining two 
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non-Gaussian states does not necessarily lead to an increase of the overall nonG. 
As 2-PHAV is a state described by 2 parameters, namely the mean value, Mr, 
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Fig. 3. Upper left panel: nonG measures ca and cb as functions of the mean number of detected 
photons of almost balanced 2-PHAVs (R = 1.24). Upper right panel: ca and cb as functions of the 
balancing between the two components of the 2-PHAV, at fixed mean number of detected photons 
(Mr = 4.12) of the overall state. Lower panel: ca and cb as functions of the mean number of 
detected photons of single PHAVs. Empty symbols: experimental data; Full circles: theoretical 
expectations. 



and the balancing R between its two components, we decided to better investigate 
its nonG nature as a function of one of these variables by keeping fixed the other 
one. In the upper left panel of Fig. [3] we show 6a and 6b as functions of the mean 
total energy of the 2-PHAV for a fixed choice of the balancing, namely 1.24: we 
can notice that the values of both the measures increase at increasing the mean 
number of detected photons. Moreover, in the upper right panel of the same figure 
we plot the lower bounds of the nonG measures as functions of the ratio between 
the two components at fixed mean number of detected photons, that is Mx = 4.12. 
As one may expect, it monotonically decreases at increasing the balancing. In fact, 
the most unbalanced condition reduces to the case in which there is only a single 
PHAV, whereas the most balanced one corresponds to have a balanced 2-PHAV. 
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For the sake of completeness, in the lower panel of the same figure we show the 
results we obtained for the single PHAV as a function of the mean value. Also in this 
case the experimental data, which are superimposed to the theoretical expectations, 
testify the accordance between the two measures since they increase their value at 
increasing the energy of the state. 

As in all the cases presented in Fig. [3] the behavior of the two measures is 
very similar except for the absolute values, we decided to test their monotonicity 
by following the suggestion of Ref. [17]. In the left panel of Fig. [4j we plot the 
experimental values of measure B as a function of those of measure A for the three 
cases presented in Fig. [3j It is evident that the two measures are monotone to each 
other, even if the absolute values are different. In particular, measure B is endowed 
with higher values, thus resulting more sensitive to small differences in the choice 
of the parameters. This property, together with the fact that €b is characterized by 
smaller error bars with respect to the other one, can be considered a good criterion 
to choose one measure instead of the other to quantify the nonG amount and thus 
discriminate the states under investigation for possible applications. Nevertheless, 
it is interesting to notice that the behavior of one amount with respect to the other 
is not described by a unique curve. In fact, in the case of the 2-PHAV at fixed 
ratio and variable total energy (blue dots) the curve is slightly different from those 
corresponding to the other two cases, namely the 2-PHAV at fixed total energy 
and variable ratio (red dots) or the single PHAV with variable energy (black dots). 
In order to prove that this difference does not depend on the reliability of the 
experimental data, in the right panel of Fig. [4] we present the results we obtained 
by a simulation. We plot the theoretical behavior of a single PHAV at different mean 
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Fig. 4. Left panel: cb as a function of ca f° r the three experimental cases presented in Fig. [3] 
Right panel: simulated behavior of c-q as a function of ca for different choices of the parameters 
describing PHAVs and 2-PHAVs (see the text for details). 



numbers of detected photons as black line, whereas we used colored squares + line 
to indicate the 2-PHAV at fixed ratio (R = 0.2,0.5,0.8) and variable total energy, 
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and colored dots + line to indicate the 2-PHAV at fixed total energy (M T = 2, 4, 6) 
and variable ratio. It is evident that there is not a unique curve, as already testified 
by the experimental data. Nevertheless, we want to notice that there are some 
limits in which the curves are superimposed (this happens either when the 2-PHAV 
is almost unbalanced or when it is very low populated as in both the cases it reduces 
to the case of a single PHAV) or intersect each other (such as in the case in which 
the 2-PHAV is characterized by a precise choice of total energy and ratio). 

4. Concluding remarks 

In conclusion, we have presented an experimental investigation of the non-Gaussian 
nature of the class of PHAVs by reconstructing their Wigner function and using 
two different measures, both based on quantities experimentally accessed by direct 
detection, to quantify the nonG amount. We proved the consistency of the dif- 
ferent approaches and tested the monotonicity of the two measures. Nevertheless, 
the comparison performed on diagonal states belonging to the class of PHAVs for 
different parameters settings showed that there is not a unique curve describing 
the behavior of one measure with respect to the other. In addition, we discussed 
the choice of the best measure between the two proposed. According to data, 6b 
seems to be better because it has higher absolute values and a reduced sensitivity 
to experimental errors. 
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